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I. DESCBIFTION OF A STOCHASTIC MODEL 
FOR COMBINATIONS OF RANDOM LOADS 



Many physical structures are threatened by combinations 
of loads of varying magnitudes from various sources. For 
instance, bridges, piers, dams, and buildings can experience 
loads frcm wind, snow, ice, tides, earthquakes, and so 
forth. In many instances the total load or stress experi- 
enced by a structure varies in time in an apparently random 
fashion. Certain components of the loads vary rather 
slowly; others occur more nearly as impulses, such as those 
associated with winds or earthquakes. The problem is to 
design a structure to withstand the superposition of random 
loads - from many sources with at least an approximately 
understood probability. 

The purpose of this thesis is to describe and investi- 
gate certain simple but somewhat realistic probabilistic 
loah models for use in design, and perhaps safety assessment 
of structures. In this thesis we confine attention tc the 
superposition of just two load types; sh oc k loads, and 
constant loads . For example, winds gusts, flash floods, and 
earthquakes have varying magnitudes and have relatively 
short durations in comparison t'o the times between their 
occurrences. These will be modelled as instantaneous shock 
loads. Cn the other hand, snow, ice, or water accumulation, 
or even the presence of slowly moving vehicles present loads 
that remain nearly constant in time, occasionally changing 
to new levels. These will be modeled as constant loads that 
change infrequently. Throughout this investigation it will 
be assumed that the effective stress exerted by several 
types of loads acting simultaneously can be expressed as a 
linear combination of the component loads, the load 
components being treated as stochastic processes. 
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The times between changes in magnitude of the constant 
load process are independently and identically distributed 
with distribution function H. The successive magnitudes of 
the constant load are independently and identically distrib- 
uted with distribution function F. Given the constant load 
process, the shock load process is a compound Poisson 
process. The successive shock load magnitudes are indepen- 
dently and identically distributed with distribution func- 
tion G. The conditional rate of arrival of shocks, given 
that tbe magnitude of the constant load process is x, is 
M(y.) . 

let X (t) be the magnitude of constant load process at 
tine t; and Y (t) be the magnitude of shock load process at 
time t; Z (t) =X (t) + Y ( t) is the superposition of the two loads 
at time t, and M(t)=^UP Z(s), the maximum load combination 
in [0,tj, see figure 1.1. 

Let Tx = {inf t>0: Z(t)>x } ie, the first-passage time 

for the load combination to exceed a stress level x. T, 
will represent the time to failure of a structure whose 
strength is x, and is subjected to a stress history 
{Z(t) ; t >03 . 

Gaver and Jacobs (1981) studied the above model for the 
case in which the distribution H is exponential and the 
shock load rate m is a constant independent of the shock 
load process. Related load combination models that have 
been studied include Peir and Cornell (1973) [Ref. 1] Wen 
(19''7) [Ref. 2] Pearce and Wen (1983) [Ref. 3]. 

Asymptotic results which appear in Gaver and Jacobs 
(1984) and Jacobs (1 984) will be described in chapter 2 for 
the distribution of Tx <x> and for the tail of the 

distribution of Tx for finite x. 

In chapter 3, a simulation program for the load 
combination model will be described. 



Q 










Figure 1.1 The Load Combination Process. 



In chapter 4, a simulation study of the use of the 
asymptotic results of chapter 2 to approximate the 
distribution of Tx will be described. 
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II. ASYMPTOTIC RESDLTS FOR THE D IST RI BPT ION OF 
THE FIR ST- PASS AGE TIME Tx 



A. THE DISTRIBDTIOM CF T^ FOE LARGE X 

This section suinniarizes some of the results found in 
Gaver and Jacobs (1981) [Ref. 4] for the model in which the 
distribution H is exponential (x) and the shock load arrival 
rate m is constant independent of the constant load process. 

First, we consider the distribution function of the 
maximum load combination to occur during [0,t], M(t). 



H(t)=P{M(t)<x3 

=P [M (t) <x, T( >t} +P [M (t) <x,T, <t} (2-1) 

where Ti is the arrival time of the first shock. 

A renewal argument yields; 

F {M (t) <x,T, >t}= 

? {M (t)<x,T, <t} = J^^Ae'^''d'/j/€xp(-x/u&tx-^j)HxCt-W/:tc(^; 
where G(x)=1-G(x). 

Next take Laplace transforms with respect to t of (2-1) . 



ly^) = io di 

=M^(g )+ AMx(? )Hx(I) (2-2) 

where 



^(? ) = //C§-^A 4 F(dy) 

M,(%) 



X 

Hence 
) = 



(2-3) 



(2-4) 



•x'^' , 

It seems to be difficult tc invert H^(f) for any inter- 
esting choice of the distributions F (x) and G (x) . However, 
useful information can still be gleaned from equation (2-4). 
First, from the definition of Tx , 

E {M (t) <x}=P (Tx >t) (2-5) 

Thus, H^(f) = p (Tx >t) dt (2-6) 

let % -*o, in equation (2-6) to find that 



m{x)=E(Tx)=Hx(0) = 



Mx(.o) 



(2-7) 



The next result is the limiting property for the first 
time the load combination process exceeds a given level x. 
The limiting distribution of Tx is exponential in the sense 
that 

liff E{m (x)'* ?x>t} =e-'‘' (2-8) 

A ^ Ao 

where x^ =inf {t ; F*G (t) = 1) 

The prccf of this result for a more general model that 
includes that of chapter 2 can be found in Jacobs (1^34) 
[Ref. 5]. 



B. THE TAIL OF THE DISTRIBOTION OF Tx FOR FINITE X 

In this section we will summarize results for the asymp- 
totic behavior of P(Tx>t) for large t and fixed x as t-^’OP 
for the model in which the distribution H is exponential and 
the shock load arrival rate is constant. ?lore details and 
results for a more general model that includes that of 



chapter 2 can be found in Jacobs (1984) . 

Frcm the equations (2-1) and (2-5) 
p (Tx >t) = f 

* tLs j^puv pcTx>t-i; (2-9) 

let 

L (ds)= // pcd^) ds (2-10) 

I(t)= (2-11, 



The renewal equation (2-9) gives; 

P (Tx>t)=g(t) + Hx(t-s) (2-13) 

where g (t) = (2-14) 

Following the development in Feller (1971; page 376) 

[Ref. 6 1 it will be assumed that there exists a number^ 
such that 

1= J* (2-15) 
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This root is unique and since the distribution of L is 
defective, is positive. 

Let 

L*(ds)= L(ds) (2-16) 

g*(t)= get; (2-17) 

(T,^>t) = (t) (2-18) 

Then 

P*(Tx>t)=g*(t)+ (t-s)}L* (ds) (2-19) 

It is assumed that the assumption for the key renewal 
theorem holds. Application of the key renewal theorem 
yields 

Lira e^"'^^P(Tx >t) = -p/"g^{t)dt (2-20) 

( t) dt= F ( dy) 

(=- 21 ) 

Rewriting the equation (2-15), the equation determining 
X s 

1= /; A ds 

= Ppedy; — ^ (2-22) 

A 4 xAs-u-yj -Xc^; 

and 

= Ifpc<^V - )icx)y (2-23) 

The key renewal theorem implies; 

Lira ^ (Tx >t) = — (2-24) 

a/,Vco(^;- ^ - 

Examples: C + ^ 6U-i; - 

It appears to be difficult to solve for j< analytically 

in general. We will discuss the numerical computation- of X 

for two examples. In both examples, and the 

distribution of the shock and constant load magnitudes F and 

G are of the form; 

F (X) = 

G (X) = 
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In example 1, a=h=1 and in example 2, a/b=2. 

To compute the j( for T(x)= e"®’^ ,^(x)= e'*’^ , 

let 

^ iX) - pcd'i) X + xt frCx-^) -;kw; (2-25) 

where A = -d =1. 

Icr model 1, 

f (x) = -j7i^{(/'X)0-e-^;-e'’'CMi'5c+e->';-i^((/-x)€-*4e'^B] — (2-26) 

for>i<1 + e"^ andxti. 

7orX=1, f (x) has a singular point, so by using the 
I'HOPITAL 3 UL^ [Ref. 7] 

^im f(x)= Lira g'(x)/h'{x) 

where g (jc) is the numerator and h(x) is the denominator 



of t]ie function f (J^) . 

For X =1 , 

f (1)=f (e^-e-^) (2-27) 

For model 2, 



+ (2-28) 



for X <1+e"‘”' and X t1. 

ForX=1 f (X) has also singular point, so by using the same 
iT'ethcd as used in model 1 ; 
f orX =1 



1 ( 1 )= 

Differentiation of the equation (2-25) gives 

f'(x)=J: 



(2-29) 






(2-30) 



Equation (2-30) is greater than 0 for 0<X< 1 +6"^^ and x . So 
the function f (x) is a monotone increasing function in the 
interval 0<X.<1+e"‘”^ . Using this result, we compute the dc 

from the equation (2-22) by applying the bisection method 
[Ref. 8]. 

Using the computed X we compute the constant such that; 

(f* = pCol^) ^ 

is the constant of equation (2-24). 



ir cr case a=b= 1 , 

-[t-^ iz-x) - ze-‘ l~ie-“iz-li>) - ]j ■■ ( 

For case a/b=2. 



1-31) 



2-32) 



III. COMPOTER SINDLATION 



This section describes the computer simulation model. 
fie denote the constant load arrival rate by X, the shock 
load arrival rate by^, the distribution of constant load 
magnitude by F (x) , and the distribution of shock load magni- 
tude by G (x) . The simulation model will be used to study 
the quality of the approximations of the distribution of Tx 
by the asymptotic results described in chapter 2. 

The following describes the simulation algorithm. The 
stress levels are x, <x^<X 3 

1, Set T=0, No=0- 

2, Set V=0, generate a constant load magnitude, C, and 
duration interval Tc. Find the largest index i>No 
such that X; <C . 

Put N, equal to that index and 
set 

Tx: =T for N, <i<K, , 
set No =N, . 

3, Generate a shock load magnitude, S, and a time until 
the shock occurs, Ts. 

4, If the shock interval does not exceed the excess life 
of the constant load duration, Tc, find the largest 
index i>No such that X; <C+S. 

Let N, be that index. 

Set Txi =T+V + Ts for No <i<N, . 

Set No =N, , V=V+Ts, Tc=Tc-Ts. Go to step 3. 

5, If the shock interval exceeds the excess life time 
of the constant load duration, 

set T=T+Tc, go to step 2. 
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One realization of the simulation is completed when one 
Tx is computed for each level X; ,1<i<n. Each simulation 
consisted of 5000 replications. Further details are 




Figure 3.1 Generating T. 
included in chapter 5. 

The hey subprograms are following. Other programs will 
be explained by reading from the output results. 

EAKDC ; generate the Tx . 

CCFPUT ; compute the P(Tx = 0), true P(Tx=0), sample mean, 
true mean, standard deviation, coefficient of 
variance. 

Sample probability 

P(Tx=0)= . T* =-g is obtained from 

the sample. 

True probability 

P (Tx=0) = 1-P(Tx>0) =1-P (M(0) <x) =1- j*F (dy) =F (t) is 
obtained from the equation (2-5) . 

The sample mean 
X =;^L Txt and 
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the true mean E(Tx,) is obtained in special 
cases from the equation (2-7) . 

See details in Gaver and Jacobs (1981). 

Sample standard deviation 
(^ = obtained from the sample, 

and 



Sample coefficient of variance 

= = is obtained from the sample 

mean and sample standard deviation. 

COM? ; compute tie sample quantiles, and the quantiles 
of an exponential distribution with mean E(Tx,). 

A sample quantile is 

= (o<'< sample size) order statistic of sample and 

id = I ~ is the quantile of exponential 

distribution with mean E (Tx ) . 

C CM ; compute th.e quantile of the sample conditional 

distribution given Tx>0 and the exponential distrib 
ution with mean E (Tx ) /P (Ty>0) . 

The sample quantile of the conditional distribut 



SK 



INT 



ion 

eC y (number of positive Tx ) order statistic 
is obtained for sample of positive Tx . 

The conditional quantile of the exponential dist 
ribution with mean E (Tx ) /P (Tx>0) is 

; compute the dC such that; 

'K={'K(x): 
bisection method. 

compute the constant such that ; 



A 

x+K 



} by applying the 



A 



C A+ K 

by substituting the 'K found in SK. 
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Ihe following cases were simulated. 



distribution 


case 


A 


MCO 


F (X) 


G (X) 


a , exponential 


1, 


1 


1 


l-e'* 


1-e'"' 


(a=b=l) 


2, 


1 




11 


11 




3, 


2 


1 


11 


11 


b , exponential 


1, 


1 


1 


l-e-^'' 


l-e'""" 


(a=2,b=1) 


2, 


1 


2 


11 


11 




3, 


2 


1 


11 


11 


c, exponential 


1, 


1 


1/c 




l-e""^ 


with varying 


2, 


1 


e"*' 


1? 


11 


arrival rate 


3, 


1 


e'^ 


11 


11 


d,Fareto case 


1, 


1 


2 




1-e'^ 




2, 


1 


2 


-j ' 


l-e'’^ 




•3 

f 


1 


2 




1-Tk- 


These different : 


Stitt ings 


will 


be referred to 


as cas 


(A-2), (A-3), (3- 


-1)/ (B- 


2) , 


(B-3) , 


(C-1) , 


(C-2) , 



(A-1) , 
(C- 3) , 

(D-1), (r-2) , (D-1) respectively throughout this paper. For 

each above case, we compare the asymptotic result (2-8) and 
the experimental results. 

The next setting is to examine the exponential approxi- 
mation of the tail cf the distribution of Tx, for finite x. 
The cases of identical exponential distributions, 

F(x)=G'(x)= e’’' and different exponential distributions, 

F (x) =e'‘*'’‘ G (x) =e'^^ where a/b = 2 
were considered. 

The following conditions are considered. 

e, F (X) =G (X) =e->^ 

f, X=M=^ F(x)=e'^'' ,G(x)=e'^ 
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These different settings will be referred to as case (S) and 
case (F) respectively. For each case, compute the , and 
constant and compare the quantiles of the asymptotic 
results (2-24) to the simulated data for each level of x. 
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IV. ANALYSIS OF THE SIMOLATION RESOLTS 

As indicated in chapter 2 the exact distribution of 
is difficult to obtain in general. In this chapter simula- 
tion will be used to study the accuracy of the approxima- 
tions of the distribution of T,(. by the exponential 
distributions suggested by the asymptotic results of chapter 

•7 



A. THE DISTRIBOTION CF Tx FOE INCREASING X 

The limiting result (2-8) indicated that the distribu- 
tion of Tx approaches the exponential as level x increases. 
This result suggests that the distribution of Tx can be 
approximated by an exponential distribution at least for 
large x. In this section this exponential approximation will 
be studied via a simulation. 

"I • Constant and Sh ock Load Ma gnit udes Iden tica l 

Exponent ia 1 Distri bu tion s 



In the first collection of three models to be 
considered F (x) =G (x) = The times between constant load 

changes are exponential with parameter A and the shock load 
arrival rate is a constant independent of the constant 
load process. 

For these cases it is possible to derive an analyt- 
ical expression for E(Tx). 

■p frr ^ p A < / - ~ ^ Jon C X + f ii — 1 1 

a^ l 1 1“ 

In the Tables, P(Ty =0) refers to the simulated 
sample probability that Tx =0; X-3AR gives the simulated 
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sample mean of ; ST-DEV is the simulated standard 

deviation of ; VAE (X-B) is the variance of the sample 
mean; COEF-V is the sample standard deviation divided by the 
sample mean. T P(?x=0) is the theoretical probability that 
Tx. = 0, F (X) ; T x-bar is the theoretical mean of Tx computed 
from the equation (4-1) . 

The simulated coefficient of variation in all three 
cases decreases as the level x increases becoming quite 
close to the theoretical exponential value of 1 when x=5. 

To compare the simulated distribution of Tx to the 
approximating exponential distribution, quantiles from the 
simulated data V7ore computed and compared to the corre- 
sponding quantiles cf an exponential distribution having 
theoretical mean (4-1) . The simulated quantiles were 
computed as described in chaj^ter 3. The exponential 
o(-quantile is given by 
q^=-a.ln (W) 

where a is the mean of the exponential. 

For each level of x, the first row in the tables 
(A . 1 ) , (A. 2), {A. 3) for the quantiles of the distribution of 

Tx gives the quantiles of the simulated data and the second 
row gives the corresponding exponential quantiles for an 
exponential distribution having theoretical mean (4-1). 

Cne way the distribution of Tx differs from an expo- 
nential is that it has an atom at 0. In particular, 

P (Tx=0) =F (X) =e'^ 
for the case considered here. 

Quantiles were used to compare the simulated condi- 
tional distribution of Ty given Ty >0 to an exponential 
distribution having mean E (Tx ) /P (Tx >0) . The simulated quan- 
tiles for the conditional distribution were computed as 
described in chapter 3. These quantiles were compared to 
those of an exponential distribution having theoretical mean 
Z(Tx) /P (Tx>0) . 
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In the tables (A.1), (A. 2) , and (A. 3) for the guan- 

of the conditional distribution for each level x, the 
row shows the simulated quantile and the second row 
orresponding approximating exponential quantile. 

The quantiles of the simulated conditional distribu- 
are much closer to their exponential approximation than 
uantiles for the unconditional distribution to their 
ential approximation. However the quantiles for the 
ditional distribution get closer to those of their 
ximating exponential as x increases. 

Comparing tables (A.I), (A. 2) , and (A. 3) , it appears 

increasing the arrival rate of shock or the rate of 
e of the constant load magnitude decreases the guan- 
. The change of arrival rate of shocks appears to have 
reater effect. 

• Co nst ant and S hock Load Ma gn itud es Have Diff eren t 
Di str ibutions 



In the next three cases, the distribution of the 
load magnitudes is exponential with parameter 1 and 
istribution of the constant load magnitudes is exponen- 
with mean 0.5. The other model assumptions in tables 
, (B.2), (B.3) correspond to those in tables (A. 1) , 

, and (A. 3) . 

As before, for each level x, the first row in the 
ile table gives the simulated quantile and the second 
ives the approximating exponential quantile using the 
ential distribution with theoretical mean. Comparing 
uantiles of T< in table (B. 1) with (A. 1) (respectively 
with (A. 2) and (B.3) with (A. 3) ) , it is seen that 
asing the mean constant load magnitudes has increased 
uantiles. Further, it appears that the convergence of 
istribution of Tx to an exponential is faster in the 
of smaller mean constant load magnitude. 



23 



The exponential approximation to the conditional 
distribution of T_x gi^ven T^(,>0 has theoretical mean 
E(Tx)/P(Tx, >0). Once again for small levels x,it appears 
that the exponential approximation to the conditional 
distribution of T* is better than that for the unconditional 
one. 

3 • Arri val Rate of Shock L oa ds Dep ends on the Con st an t 
ioii Ma gni tu de 



In the next 3 cases, constant load magnitudes and 
shock load magnitudes are exponential with mean 1. Constant 
load magnitudes change according to a Poisson process with 
rate 1. Given the constant load magnitude at time t is C, 
the probability a shock load v.’ill arrive in the time 
interval [t,t+h] is x^(C)h+o(h) . 

In case (C-1) the conditional shock load arrival 
rate is/U(c)'*=C. Comparison with table (A.1) indicates that 
conditional shock arrival rate x(c)~* =C has the effect of 
decreasing the quantiles of Tx . Further the exponential 
approximation to the quantiles of Tx is not as good as in 
case (A-1) . 

The approximating exponential distribution to the 
conditional distribution of Tx given Tk>0 has a mean of 
F (Tx ) /P (Tx >0) . The quantiles of the exponential approxima- 
tion to the quantiles of the conditional distribution appear 
to be closer than those for the unconditional distribution. 

In case (C-2) shock loads arrive with conditional 
arrival rate ^C)"* =exp(C). Comparing the table (C.1) and 
(C.2), it is seen that the quantiles of the case ^(C)"^=C are 
less than those for the case M.(Cf'=exp (C) . This is because 
interarrival times tends to be larger for the case 

=exp(C). Comparing tables (A.1) and (C.2) indicates 
that the exponential approximation to the quantiles of ""x is 
not as good as in the case /H=1. 
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In case (C-3) the conditional shock arrival rate is 
=exp (-C) . Comparing tables (C. 2) and (C.3) indicates 
that the guantiles fcr the case 

yU.{C)'^ =exp (-C) are less than those for the case =exp (C) 

and furthermore the exponential approximation of Tx appears 
to be close for the case (C-3) than that of for the case 
(C-2). In all of the cases the simulated mean was used as 
the parameter in the exponential approximations. 

^ • 10^3 Ma^ni^des Rave a Pareto Di stri but ion 

In the next three cases, the times between constant 
load changes of magnitude are exponential with mean 1. Fhock 
loads arrive according to a Poisson process with rate 2. In 
case (D-1) the constant load magnitudes have a Pareto 
distribution with parameter °<=1 and the shock load magni- 
tudes are exponentially distributed with mean 1. Comparison 
with table (A. 2) indicates that the Pareto constant load 
magnitude has the effect of decreasing the juantiles of Tx- 
Further the exponential approximation to the guantiles of Tx 
is not good in the Pareto case; (the approximating exponen- 
tial has a mean of the simulated sample mean). 

The approximating exponential distribution of Tx 
given conditional distribution has a mean of 

F (Tx ) /P (Tx >0). Once again the guantiles of exponential 
approximation to the guantiles of the conditional distribu- 
tion appear to be closer than those for the unconditional 
distr ibution . 

In case (D-2) the constant load magnitudes have a 
Pareto distribution with parameter c< =2. Comparing tables 
(D.1) and (D. 2) it is seen that the guantiles of the simu- 
lated distribution of Tx for the caseot=2 are larger than 
those for the caseo<=1. Comparing tables (D.2) and (A. 2) 
indicate that the exponential approximation to the guantiles 
of Tx is not as good for the Pareto case as for the 
exponential case. 
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In case (D-3) the constant load magnitudes have 
Pareto distribution with parameter « =1 and the shock load 
magnitudes have Pareto distribution with (X=1. In all three 
cases the simulated mean was used in the approximating guan- 
tiles . 



B- THE EXPONENTIAL TAIL OF THE DISTRIBUTION OF T^ FOR 
FINITE X 

In this section simulaticn will be used to study the 

exponential approximation 

suggested by the asymptotic result 

/ 

^ • 



lini =-^ 



This is an approximation for P (Ty >t) for fixed finite x; it 
sliould become more accurate as t becomes large. 

Two cases are considered. In both cases shock loads 
arrive according to a Poisson process with rate 1 and 
constant load magnitudes change at the times of arrival of a 
Poisson process with rate 1. The shock load magnitudes have 
exponential distribution with mean 1. In case (E) , the 
constant load magnitude distribution is exponential with 
mean 1; In case (F) it is exponential with mean 0.5. 

Description of the computation of X und for these two 
cases can be found in chapters 2,3. It can be seen from the 
computed value of X(x) and ^(x) appearing in tables, (E.1) 
and (F.1) that JC(x) is approaching the value 1/E(Tx) (where 
E (Tv ) is the theoretical mean) as x becomes larger and for 



Further , y (x) is approaching 



the cases computed X(x)<t^, 

1 as X becomes larger. 

. To study the exponential approximation to the distribu- 
tion of T , quantiles from the simulated data were computed. 
These can be found in tables (E. 2) and (F. 2) . For each level 
X, the first row presents the simulated quantile; the second 
row gives the approximating c<-guantile computed by 

I r 






( I ) 
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The third row shows the approximating -quantile computed by 

qo(=-E ) In (1-c<) ( n ) 

where the theoretical mean is used for E (T^ ) . 

Comparing the quantiles in tables (E. 2) and (F.2) it can 
be seen that the approximating quantiles (X) are close to 
the simulated ones even for The approximating quan- 
tiles computed in (i) do less well but improve as x becomes 
larger. The two approximating quantiles approach one another 
as X becomes larger, as expected. Approximating quantile (]L) 
has the advantage, ho wever, of being easier to compute. 

Table (G.1) presents 90 T confidence intervals for 
o 4 -guantiles with <>(.>0.9 for case (E) . These confidence inter- 
vals were computed using a large sample approximation, see 
Conover (1980; page 111-114) [ Bef . 9]. 

Comparison with table (E.2) suggests that the approximating 
quantiles approximates the simulated ones well for these 
values. 



C. CCNCIDSIONS 

Approximating quantile X is always better than approxi- 
mating quantile £ . Approximation i approximates q^ well for 
0 ^. as small as 0.40. However i is more difficult to compute 
than I. . 



The performance of approximating quantile IL can be 
improved by chancing it to 

PCTx>o) 

and using it to approximate the quantiles of the conditional 
distribution of Ty given Ty>0; for the models considered 
P{"')(=0) is easy to compute being equal to F (x) . 
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TABLE A . 2 

Identical Exponential( case 
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V. SIHOI^TION A LGORIT HM 

A. VARIABLE DEFINITICH 

Tx; ={inf t>0 ; Z(t)>Xi } 

where Z (t ) =superposition of constant load and shock load 

Nrow = iraximum number of level X; ; X,<X^<X^< 

Ncol=maximum number of repetition of generating Tx. 

Lc=constant load magnitude 

Ls=shcck load magnitude 

Tc=constant load interval 

Ts=shcck load interval 

Ipc=ccunting the number of repetition 

I = counting the maximum level of X; 

E. ALGORITHM 

Input ;parameters { a, b, seed numbers ) 

Initialize the stress level X;; X, <X^<X^< 

i pc=0 
Repeat 

set Tc=0 
M=0 
TTs=0 
Repeat 

generate constant load;Lc 
generate constant interval;Tc 
Rind L such that L= {sup i;Lc>x; } 

Set Tx: = TO for 1<i<L 
M=L 
Ci=Tc 

A, generate the shock load;Ls 
generate the shock interval ;Ts 

B, If Ts>Ci, then 



45 



TO=TO +Ic 
Ts=Ts-Ci 
TTs=0 

generate constant load;Lc 
generate constant interval; Tc 
Find 1 such that L= {sup i; Lc>Xi} 

If L>M, then 

Set Tx; =T for M<i<L 
M=L 

Else continue 
Ci=Tc 
Go to E 

Else, Find 1 such that L={sup i;Ls+Lc>Xi } 
TTs= Ils+Ts 

Set Tx; =Te+TTs for H<i<L 
Ci=Ci-Ts 
Go to A 
Until (L=Nrow) 

Ipc=Ipc+ 1 
Until (Ipc=Ncol) 

End Algorithm 
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WRITE (fc , 125) ( X( I) ,I = 1,NR0W) 

125 FORMAT (/, IX, • (X), I/EITX), (X) 
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